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Abstract 

Far tails of the density of state (DOS) are calculated for the simple models of organic amorphous material, 
the model of dipolar glass and model of quadrupolar glass. It was found that in both models far tails are 
non-Gaussian. In the dipolar glass model the DOS is symmetric around zero energy, while for the model of 
quadrupolar glass the DOS is generally asymmetric and its asymmetry is directly related to the particular 
geometry of quadrupoles. Far tails of the DOS are relevant for the quasi-equilibrium transport of the charge 
carriers at low temperature. Asymmetry of DOS in quadrupolar glasses means a principal inequivalence of 
the random energy landscape for the transport of electrons and holes. Possible effect of the non-Gaussian 
shape of the far tails of the DOS on the temperature dependence of carrier drift mobility is discussed. 
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I. INTRODUCTION 


One of the fundamental characteristics defining charge transport properties of any disordered 
material is a density of states (DOS) P(U), i.e. a distribution density of random energy U of a 
charge carrier hopping in the bulk of the material. In some situations the very functional form 
of P(U) could be used for a quick estimation of the relevant transport properties of the mate¬ 
rial. For example, if the DOS decays fast for U —y — °°, then the non-dispersive quasi-equilibrium 
transport regime eventually takes place for t —» °° and the stationary density of occupied states 
P(f/)exp( — U / kT ) is developed (we assume here the Boltzmann statistics). Let the maximum 
of this distribution is located at U eq . In quasi-equilibrium regime the limiting step for the charge 
transport is a carrier escape from deep states where carrier hops up in energy to perform a corre¬ 
sponding transition in space. It is reasonable to assume that the final energy of such hop is located 
somewhere near the maximum Uq of the DOS. Hence, a rough estimation of the temperature de¬ 
pendence of the carrier drift mobility fi{T) for low electric field should be 


/i exp 


/ Uo-Ueq \ 


( 1 ) 


with Uo — U eq serving as an effective activation energy. For the Gaussian DOS 

P(l ' )= (2^) eXP ("^ 1 ) <2) 
typical for amorphous organic semiconductors [1-4], U eq — —o 2 /kT and fl °< exp [— a 2 / (kT) 2 ]. 
This simple estimation gives an exact leading asymptotics for the low field mobility temperature 
dependence in ID case [5] and differs from the corresponding asymptotics in higher dimensions 
by the numeric factor ~ 1 in the exponent [1, 6, 7]. More refined approach invoking the conception 
of so-called transport energy provides even better description of the dependence jJ.(T) [8, 9], at 
least for the materials with spatially non-correlated random energy landscape. 

Experimental data on the mobility temperature dependence suggest that in organic materials 
a ~ 0.1 eV [1] and, hence, at the room temperature o/kT ~ 4 — 5. For some materials charge 
transport has been observed in time-of-flight experiments even for o/kT ~ 6 — 7 [10], though in 
that case it is highly dispersive. This means that typically the maximum of the occupied DOS, 
relevant for the quasi-equilibrium transport, is located at the tail of the initial DOS. 

There is a strong evidence that in organic materials a major contribution to the total DOS has 
the electrostatic origin: it is produced by the interaction of charge carrier with randomly located 
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and oriented static dipoles and quadrupoles. Estimation of the electrostatic o taking into account 
typical dipole and qudrupole moments and concentration of polar molecules provides values in 
reasonable agreement with experimental data [7]. Models of dipolar glass (DG) and quadrupolar 
glass (QG) have been suggested for description of random energy landscape in amorphous or¬ 
ganic materials. In these models it is assumed that there is no correlation between orientations 
of dipoles or quadrupoles. The models naturally produce highly spatially correlated random en¬ 
ergy landscape necessary for realization of the specific Poole-Frenkel mobility field dependence 
ln/i E 1 / 2 [5, 11, 12]. 

Until now an accurate calculation of the behavior of the far tails of the DOS in DG and QG 
models has not been carried out. For DG and QG models the central peak of DOS has a Gaussian 
shape if concentration of dipoles and quadrupoles is not too low [12-14]. This shape is essentially 
guaranteed by the Central Fimit Theorem which is not valid for the far tails of the distribution. 
Typical distribution of dipoles relevant for the main body of DOS and for the far tail are qualita¬ 
tively different (see Fig 1 (a) and (b), correspondingly). In this paper we are going to carry out a 
direct calculation of the functional form of the tails of DOS in DG and QG models and estimate a 
possible effect of the tail shape on charge transport. 


II. DENSITY OF STATES IN THE DG MODEL 


Let us consider the simple cubic lattice with sites occupied by randomly oriented dipoles having 
dipole moment p. We assume that there is no orientational correlation between dipoles, their 
orientations are uniformly distributed, and fraction of sites occupied by dipoles is c. Then the 
distribution density of the random energy U is 


oo 

P(U) = U(u-Y J uX\ = ^jdyL*viy(u-Y J U, 


1 

4 7T 


dyexp(iyU)Yl 


c d (cos #„) exp(— iyU„) + 1 — c 


-l 


( 3 ) 


( 4 ) 


where angular brackets mean an average over positions and orientations of dipoles, index n runs 
over all lattice sites apart from the initial reference site and U„ is the contribution of the point 
dipole at the site n to the total energy U of the charge carrier at the reference site 


ep n r n ep 

U n =- ZTT- = -—^COSldn, 


erl 


en 


( 5 ) 
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FIG. 1. a) Typical distribution of dipoles that generates random energy 1/ in the main body of the DOS 
(random energy is measured at the point in the center of the picture), b) Large stochastic cluster of dipoles 
relevant for the far tail of the DOS. 


here £ is a dielectric constant and a is a lattice constant [14]. After integration over in equation 
( 3 ) we obtain 

P(U) = ±-[ dy exp [iyU + S(y)], S(y) = £ln ( c ^l + \ - C Y z n = ^f j- (6) 

Z7T ,J oo n \ Zn J ££„ 
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Let us calculate the integral (6) using the saddle point method. The structure of the exponent in 
the integral (6) indicates that the saddle point is located at y s — ioo where ft) is a real number. We 
will see that for far tail U —ft) —too. We obtain a rough estimation of ,S(ft)) replacing the 
summation in equation (6) with the integration 


S((o) ~ 2 KzJ 2 j ^2 In 
o 


3/2 f dz,J c sMn +l _ c y Za = (0Ud = i ^ 


For ft) —>• oo 


5(a)) = InzJ 2 


A d~ 


dz , ( sinhz . 

■In ( c -h 1 — c 


75/2 


Za 


2nz 3 J 2 A d -4nz a , 


here we keep the leading correction only and 


(7) 


( 8 ) 


A d — 


dz , ( sinhz i . 
In ( c ——-1-1 — c ) . 


7 5/2 


(9) 


We find the position of the saddle point from the equation (note that ft) > 0 is appropriate for 
U> 0) 

dS I/O / IT dir. fTj \ 2 

( 10 ) 


- u+ ^=~ u+ - 4lrt/ "=°’ =(li®' 


and the saddle point approximation for the integral (6) is 


P(U) 


1 


2it B) 


1/2 


exp [—(oU + £(«>)]. 


(ID 


Substituting the solution of equation (10) in equation (11), we obtain the far tail asymptotics 


P(U) 


t/| +4jrt/<A 

) 


1/2 


exp 


(\U\+47cU d y 


\U\ — > oo, 


( 12 ) 


in this form the asymptotics is valid for both positive and negative U. We see that far tails are 
non-Gaussian but the DOS is still symmetric around U — 0. 

Saddle point approximation is valid if the correction to the quadratic term in the expansion of 




- 1/2 


\dco 2 J 


. Estimating the correction 


S(ft>) around the maximum is negligible at the scale 8(0 ~ 
using the next term of the expansion, we obtain the necessary condition for the validity of the 
approximation as 


d 3 S 

da 3 


°c z a 3,/4 -C 1. This inequality is valid if \U\/U d 1. In fact, 


all other corrections are negligible, too. 
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III. DENSITY OF STATES IN THE QG MODEL 


In the QG model we replace dipoles with quadrupoles. A point quadrupole generates the elec¬ 
trostatic potential 



(13) 


where Q is a symmetric traceless tensor £ Q n — 0. It may be transformed to the diagonal form 

i 

Q = diag(—<2(1 + ce)/2,—<2(1 — a)/2, <2) having two scalar parameters Q and a. Any quadrupole 
could be considered as a linear combination of axial and planar quadrupoles Q = Q fl — aQ p /2, 
where Q a = diag(—<2/2, —Q/2,Q) and Q p = diag(<2, — <2,0). Let us consider two separate cases, 
i.e. the cases of pure axial (AQ) and planar (PQ) quadrupoles. 

Note that there is a very distinct difference between density of states in AQ and PQ glasses. For 
PQ the inversion Q—t —Q that transforms the positive carrier energy to negative one and vice versa 
is equivalent to the spatial rotation of the quadrupole. All spatial configurations of the quadrupole 
have the equal weight, hence, the DOS is symmetric around U = 0. For AQ the inversion is 
not equivalent to any possible rotation and the DOS is not symmetric. Substantial asymmetry is 
developing only for far tails of the DOS because for typical high concentration of quadrupoles the 
main body of the DOS has a Gaussian form and, hence, is symmetric. Asymmetry of the DOS for 
any arbitrary QG (being, in general, a mixture of AQ and PQ) reflects the particular geometry of 
a quadrupole and is directly related to the contribution of the axial component. This is not so for 
DGs, where the inversion of the dipole moment is equivalent to the spatial rotation and the DOS 
is exactly symmetric. 

A. Axial quadrupoles 

DOS can be calculated using equation (3) with the dipolar energy replaced by the quadrupolar 
one 



(14) 


and 



n 


X — COS ■&„. 


(15) 
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Again, at the saddle point y s = i(d. Replacing summation over n with integration over z — we 
obtain 


Za 

S{a >) ~ ^ z a J 3 In [cF a {z) + 1 
o 



Za — | G) | U a 


eQ\a>\ 
4£a 3 ’ 


l 

F a (z) = j dxex p [-z{3x 2 - 1)]. 

o 

For |f/| —>■ °° saddle point is located at |ft)| —> Let us consider the case (0 —> Then 


(16) 

(17) 


1 />r\ 1/2 

^ <i8) 

Hence, keeping the major contribution and leading correction only 

4jr 

S{co) « — Za (in Za +Aq) , (19) 


where 


dz 


dz . 


A q = / -j-In [cF a (z) + 1 —■ c] + / 3 ln[cF a (z)e ^(l-c)? : . 


J z 


Saddle point equation is 


4 n 


U + — (lnz fl + l) t/,a — 0 


( 20 ) 


( 21 ) 


(we see that U > 0 for (O > 0), and the final asymptotics is 



Analogous calculation for co —» —°o provides the asymptotics for f/ —> — °°. Here the corre¬ 
sponding limit for the function F a {z) is 


Fa(z) -A <? 2z , z = cou m 
6 z 

4-7T 

S{(d)tt —Za (2In Za + Bq) , 


and finally 


P(£7) ~ 


Pzjy 2 

4irt/„ 
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ex P 1- ^-Za I , Za = exp I - 


3 U 
87r U„ 


— B q /2 — 1 ) , U^~c 


where 


B q = J * In [cF fl (-z) + 1 - c] + y ^ In [cF a (-z)e + ( 1 - c)e 2z ] . 


(23) 

(24) 

(25) 

(26) 


o 1 

We see that for axial quadrupoles the far tails of DOS are asymmetric in agreement with the general 

consideration in the previous section. 
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B. Planar quadrupoles 


For planar quadrupoles 


U n — u n (cos 2 id n — sin 2 cos 2 <p n ), 


U n — 


eQ 

2 er3’ 


and 


S(y) = L ln 


1 2k 


c 

2k 


j dx j dcp exp [iyu n (x 2 (1 + cos 2 cp) — cos 2 <p) ] + 1 — c 


o o 

For y s — i(0 and a> —* °° the corresponding function (z = (0u n ) 

1 2k 


F p (z ) = — j dx I d(p exp [—z (x 2 (1 + cos 2 <p) — cos 2 (p ) ] 


o o 


2k 


4 (zk) 


1/2 


d(p 


exp(zcos 2 <p) e z 
(l+cos 2 ^) 1 ^ 2y/2z 


For large oo the asymptotics of S(co) is 


S{Za) ~ ^ (in Za "FC^) , Za — 0)Up — 


(27) 


(28) 


(29) 


(30) 


(31) 


(32) 


2ea 3 ’ 

where 

1 oo 

Cq = J ^ In [cFp(z.) + 1 - c] + j ^ In [cF p (z)e _z +(1-c)e _z ] 
o ^ t 

and this asymptotics is valid for <u —>■ ±o° if we define z a = \(o\U p in full agreement with the 

general symmetry of the DOS. Hence, the asymptotics of P{U) for planar quadrupoles is described 
by equation (22) where U is replaced by \U\,U a by U p , and A q by C q . For planar quadrupoles this 
asymptotics is valid for both positive and negative U. 


Analogous to the DG model, saddle point approximation is valid for axial and planar quadrupoles 


if 


d 3 S 


dad 


(j^|) ^ !/ 2 « 1. This inequality is again valid if \U\/U a , p » 1. 


Arbitrary quadrupole is a linear combination of axial and planar quadrupoles. For the cor¬ 
responding QG model the actual far tail asymptotics is determined by the contribution of the 
component providing the slowest decay of DOS. 


IV. DISCUSSION: IMPLICATION FOR CHARGE CARRIER TRANSPORT 

Main body of DOS in DG and QG models is Gaussian and, thus, symmetric. This means that 
the effect of the energetic disorder on the transport of electrons and holes is exactly the same if 
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the transport is determined by the main body of DOS. Analytic results and computer simulation 
show that the tails of DOS at \U\/o ~ 4 — 5 are still Gaussian [13, 14]. Far tails in the QG 
model becomes asymmetric and provide inequivalent environment for carriers of the opposite 
signs. Unfortunately, experimental observation of the quasi-equilibrium charge transport for low 
temperature o/kT ~ 7 — 8 is extremely difficult. 

Still, we could provide some estimations for the low field mobility temperature dependence 
in the DG model for low temperature. Using the approach described in the Introduction and 
estimating the effective activation energy using the position of the maximum of the occupied DOS, 
we obtain for the leading asymptotics 



(33) 


This estimation is valid when Ud/kT 1 and the second term in equation (33) is dominating. For 
the axial QG model the corresponding relation is (for one particular side of DOS) 



(34) 


with the trivial modification for the other side of the axial QG DOS or for the case of planar QG. 
In all cases the temperature dependence becomes weaker in comparison with the Gaussian DOS 
due to the more faster decay of the DOS tail. 

In fact, weakening of the dependence j u(T) for the time-of-flight experiments at low temper¬ 
ature has been observed previously [15, 16], but the experimental photocurrent transients are so 
dispersive that this effect is better attributed to the non-equilibrium transport. Borsenberger et al. 
[15] showed that in the popular Gaussian Disorder Model transition from the non-dispersive to dis¬ 
persive transport leads to the change of the low field dependence from ln/i — {lo/'ikT) 2 
to ln/i « —(<t/2 kT) 2 , e.g. the slope of the ln/i vs 1 /T 2 dependence becomes smaller by the fac¬ 
tor ?a 2. In some situations in the low temperature region the decrease of the slope becomes much 
greater (see Fig. 11 in [17]), indicating a possibility of the deviation of the shape of the DOS from 
the Gaussian one. 

Measurement of the stationary space charge limited current [18] or technique of charge car¬ 
rier extraction by linearly increasing voltage (CELIV) [17] makes it possible to study the quasi¬ 
equilibrium transport even for very low temperature a/kT ~ 10. Unfortunately, analysis of the 
transport data is not so straightforward as in the case of time-of-flight data. 

Another important problem is the modification of the mobility field dependence for the non- 
Gaussian DOS, it will be considered in a separate paper. 
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y. CONCLUSION 


We consider the shape of the far tails of the DOS in the models of amorphous organic materials, 
suitable for the description of charge carrier transport in polar and nonpolar materials (DG or QG 
model, correspondingly). It was found that for both models the shape of the tails becomes non- 
Gaussian, but there is a principal difference between two models: in the DG model the DOS is 
exactly symmetric around U — 0 while in the QG model the DOS is generally asymmetric leading 
to the inequivalence of the energetic landscape for transport of electrons and holes. 

Deviation of the DOS shape from the Gaussian one naturally leads to the alteration of the mobil¬ 
ity temperature dependence for low temperature: the dependence p{T) deviates from the genuine 
Gaussian dependence ln/i — (o /kT ) 1 2 3 4 5 6 7 8 9 and becomes weaker (see equation (33) and (34)). Ob¬ 
servation of such dependence could be a serious argument in favor of the non-Gaussian DOS in 
the tail region. Unfortunately, any feasible transport experiment should be carried out at the ex¬ 
tremely low temperature where time-of-flight transport becomes highly dispersive even for rather 
thick transport layers. For the quasi-stationary methods such as CELIV or space charge limited 
current-voltage measurements an extraction of the carrier mobility from the experimental data is 
not so straightforward. 
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